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Abstract 

Starting from a Lagrangian, electromagnetic field in the presence of a nonlinear dielectric medium is 
quantized using path- integral techniques and correlation functions of different fields are calculated. 
The susceptibilities of the nonlinear medium are obtained and their relation to coupling functions 
are determined. Finally, the Casimir energy and force in the presence of a nonlinear medium at 
Q H ' finite temperature is calculated. 

03 1 PACS numbers: 12.20.Ds, 03.70. +k, 42.50.Nn 

£■ 

> 
in 

CO 
CN 

d 



X 



* fardin ' kh@ph ys.ui.ac.ir| 
^ amooghorban@sci . ui . ac . ir 
1 msoltani@phys.ui.ac.ir 



1 



I. INTRODUCTION 



One of the most direct manifestations of the zero-point vacuum oscillations is the Casimir 
effect. This effect in its simplest form, is the attraction force between two neutral, infinitely 
large, parallel and ideal conductors in vacuum. The effect is completely quantum mechan- 
ical and is a result of electromagnetic field quantization in the presence of some boundary 
conditions. The presence or absence of boundary conditions cause a finite change of vacuum- 
energy which its variation respect to the distance between the conductors gives the Casimir 
force [1] . The general theory of forces between parallel dielectrics was worked out by Lifshitz 
et al Q. Some years later, the theory of the Casimir effect was rederived by Schwinger 
and later on this theory extended to systems consisting of some dielectric layers j^l-jlj] and 
precision experiments on measuring the Casimir force performed 6]- 9]. 

Generally, calculations of Casimir forces are based on two different approaches 10J. In 
the first approach, the total energy of discrete quantum modes of the electromagnetic field 
is calculated and usually a regularization procedure is required, although there is still no 
consensus on the regularization procedures, which often lead to inconsistent results. On the 
other hand, this formalism despite of its simplicity, is restricted to systems where energy 
eigenvalues are known, but geometries in which energy eigenvalues are known exactly are 
few. 

The second approach is based on a Greens function method where using the 
fluctuationdissipation theorem, the electromagnetic field energy density is linked directly 
to the photonic Greens function. 

In recent years considerable attention has paid to the Casimir effect because of its wide ap- 
plications in different areas of physics such as quantum field theory, gravitation and cosmol- 
ogy, atomic physics, condensed matter, nanotechnology, and mathematical physics 
On the other hand, during the last two decades, the emergence of periodically structured 
optical materialscommonly called photonic crystals- has lead to substantial progress in the 
science and technology of optics and photonics [l5]. Photonic crystals [if]]- 13| provide new 
opportunities for enhancing and controlling nonlinear optical processes. An important step 
in this direction is to provide a quantum theory of light that takes into account nonlinear 
processes and at the same time include the absorption in the nonlinear material. Several 
attempts have been made to provide a quantized theory of macroscopic nonlinear electro- 
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dynamics [l8|-[23| and most of them are focused on strictly lossless materials. Recently, a 
consistent approach that includes absorption and dispersion has been formulated within the 
frame work of quantum brownian motion [24J]. This motivated us to investigate the quan- 
tization of electromagnetic field in the presence of a nonlinear dielectric medium and also 
consider the nonlinear effects of the medium on Casimir effect which can have applications 
in both fundamental science and engineering. 

In the 1990s, Golestanian and Kardar developed a path-integral approach to investigate 
the dynamic Casimir effect in a system which was consisted of two corrugated conducting 



plates surrounded by quantum vacuum 



25|-[26|. Emig and his colleagues also used the 



path-integral formalism to obtain normal and lateral Casimir force between two sinusoidal 
corrugated perfect conductor surfaces [^tJ-J^. Recently, this formalism is merged with a 
canonical approach to calculate the Casimir force between two perfectly conducting plates 
immersed in a magnetodielectric medium {29]. In the present work the electromagnetic 
field in the presence of a nonlinear magnetodielectric medium is quantized in the framework 
of path-integrals and as an example the mutual Casimir force between conducting parallel 
plates, immersed in a nonlinear dielectric medium, is calculated. 

The layout of the present work is as follows: In Sec. II, a Lagrangian for electromagnetic 
field in the presence of a nonlinear medium is proposed and quantization is achieved via path- 
integrals. In Sees. Ill and IV, the linear and nonlinear green's functions of the electromagnetic 
field and medium are obtained and the nth order susceptibility of the nonlinear medium 
which satisfies the kramers kronig relations is determined. In Sec.V, Casimir force in the 
presence of a linear and nonlinear medium for zero and finite temperature is calculated. 
Finally, we discuss the main results and conclude in Sec. VI. 



II. FIELD QUANTIZATION 



Quantum electrodynamics in a linear magnetodielectric medium can be accomplished by 
modeling the medium with two independent reservoirs interacting with the electromagnetic 
field. Each reservoir contains a continuum of three dimensional harmonic oscillators de- 
scribing the electric and magnetic properties of the medium 3o|]- 32]. In this section we 



follow the idea introduced in 



24| to quantize the electromagnetic field in the presence of 



a nonlinear medium and for simplicity we restrict our attention to a nonmagnetic medium 
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but the generalization to a magnetodielectric medium is straightforward. For this purpose 
let us consider the following total Lagrangian density 

£ = t-'EM + £"tnat + £-int (1) 

where Cem is electromagnetic field Lagrangian density 

1 2 B 2 

Aem = ^£oE - — . (2) 

The physical fields can be written in terms of the potentials as E = — ^ — V0 and B = V x A, 
and for simplicity we work in the axial gauge where (f> = 0. The dielectric medium is modeled 
by a vector field X(o>), which is suppose to describe its electrical properties 

f°° 1 • 1 
C mat = J du-X. 2 (u, x) - -w 2 X 2 (a;, x). (3) 

Now we define i'th component of the polarization field of the dielectric medium as 



oo roo 



Pi(x) = / duv^^X^u.x) + / du du'is (2) (uj,uj')X i (uj,x)X i (uj',x 



J J 

OO /"OO poo 



+ / dw / dw' / du"u^ ) {u,u',u")X i (u,x)X i (u',x)X i {u",x) + ... (4) 
Jo Jo Jo 

In these expression, the index i can take on the values x, y and z and 
are the coupling functions between the medium and the electromagnetic field. As it can be 
seen from (Tj0), the coupling-tensor describes the linear contribution of the interaction 
and the sequence v^ 2 \ u^ 3 \ .... describe, respectively, the first-order, second-order and higher 
orders of non-linear interactions. The coupling-tensors u^\v^\u^\ ... in (j3J) are the key 
parameters in this quantization scheme and we will see that the susceptibility functions of 
the medium can be expressed in terms of these cou pling -tensors. 



The interaction part of Lagrangian is defined by 24 ] 



Cint = A P. (5) 

Now we quantize the theory using the path-integral techniques. For this purpose let us start 
with the following generating functional 33] 



Z[J] = J V[ip]expi J d A x[C{ip{x)) + J(x)<p{x)], (6) 

where tp is a scalar field and J acts as a source or an auxiliary field and different correlation 
functions can be found by taking repeated functional derivatives with respect to the field 



J(x). The above partition function is Gaussian since the integrand is quadratic in fields. 
To obtain the generating function for the interacting fields, we first calculate the generating 
function for the free fields 

Z [Jem,{Ju}\ = J V[A]V[X] 

x exp [i J d 4 x{C EM + C rnat + J EM • A + J dujj w - X}] . 

(7) 

Using the 4- dimensional version of Gauss's theorem 

/ d 4 xC em = - J d 4 x(A-(y x y x A) + A-d 2 A), (8) 
and using the integration by parts 

/ d 4 x X(a>,x) • X(w,x) — — I d 4 x X(w, x) ■ -^-X(a>,x). (9) 
From Eqs. ([8]) and (Q, the free generating functional ([7j) can be written as 

Z [3em, Jw] = J V[A]V[X] exp - % - [ J d 4 x{(A • y x y x A + A • d 2 A) - J EM ■ A 
dujX(oj, x)(— + u 2 )X(x, u) - J u {x) ■ X(x, to)}] 

(10) 

The integral in Eq. ( jTUl) can be easily calculated from the field version of the quadratic 
integral formula and the result is 

Z [Jem, J(w)] = exp -- [ d 4 x d x'J EM ■ G EM (x - x') ■ J EM 

+ Jd 3 xjdtj dt' J dwJ u (x,t) -G u (t-tf) ■ J w (x,0], 

(11) 

here the space component of the point x G M 4 is indicated in bold by x 6 M. 3 and the 

«(o) 

time component by t or x E R. The Green functions — £') and G EM (x — x') are the 



propagators for free fields and satisfy the following equations 

d 2 ,tt(°) 

at 

a 2 



(y x y x + — )G EM (x-x') = 5{x-x') (12) 



+ u 2 }G^{t - t') = 5(t - 0<W- (13) 



Taking the Fourier transform of Eqs. (Tl2|) and (1131) we find 



G 



■MW 

EM,af3 



and 



o(0)± 

^EM,a/3( k . W ) 



5 a p - kgkf, 
k 2 -co 2 



1 



w 2 - a;' 2 - z0+ 



(14) 

(15) 
(16) 



where k is a unit vector along k and G and G refer to the longitudinal and transverse 
parts of Green's tensor in Fourier space. For further simplicity we define 

Jp = J dwu^\u)3 u + J dcu J du'v&^UjU^JuJuf + J dtu J dco'u^ 3 \ij0, u', u") 3 U J U /J u ir + • 



(17) 



Now the generating functional of the interacting fields can be written in terms of the free 



generating functional as 



34] 



n=0 ' ^ J J 

+ / dwi / du 2 ^ 2 \oJi, 002) 
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SJem(z) dzo83 Ul (z) 
8 8 8 



+ / dui / du 2 / gL; 3 z^ (3) (a;i,a;2, UJ3) 



SJem(z) dz 83 Ul {z) 83 U2 (z) 
5 8 



83 EM (z) dz 83 W1 (z) 83 W2 (z, ) 83 UA (z 



+ 



Z [J 



EM, 



(18) 



where Z[0] a normalization factor. 



III. LINEAR GREEN TENSORS 

In this section we consider a linear medium, so the sequence of nonlinear coupling functions 

^(2)^(3) 

, etc. are zero. Following the procedure presented in [29|], the linear Green's function 
of the electromagnetic field can be obtained as 



^ (l - y) = i ^(y) l^ (19) 
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To facilitate the calculations we take the time-Fourier-transform of the Green's function 
f|T9|) . after some lengthy but straightforward calculations we find the Green's function in 
frequency-space as 

G EM (x-x,u) = G (x-x,w) 

+ J ^[^(x-xx.w) / ^'{z/ 2 (o; , )w 2 G^(a;)}G'^i / (x 1 -x»]. (20) 

It can be shown that this Green's function satisfies the following equation 

V x V x G EM (x - x, u) - u 2 e(u)G EM {x. - x, u) = S(x - x') (21) 



with the following formal solution 



C S (k,-) = -J^y (22) 



G 

where e(u) = 1 + x^fa) * s the linear permittivity of the medium in frequency domain and 
the linear susceptibility x ( w ) can be written as 

,»> M = (2 4, 

Jo a/" 5 — a; 2 — z0+ 

These are complex functions of frequency which satisfy Kramers-Kronig relations and have 
the properties of response functions i.e, e(— w*) = e*(oj) and Ime{uS) > provided that 
(i/W) 2 (— u*) = (u^) 2 (u). It can be shown that these functions have no poles in the upper 
half plane and tend to unity as u — > oo. As a consequence, the electric susceptibility of 
the medium in the time domain can be written as 

Jo u 

(25) 

t < 



x (1) (*) 



which is real and causal. This proves that the constitutive relations of the medium, with 
arbitrary dispersion and absorption, is properly described by the present model. If we are 
given a definite permittivity for the medium, then we can inverse the relations (12^1) and find 
the corresponding coupling function (u) as 



u V>(u) = J-/mxW(oj). (26) 



In a similar way we can obtain the other correlation functions among different fields 29| . 
but here we consider only the following correlation function 



g em,p( x_x ' u ) =iux ( 'MG EM (x-x',w), 



(27) 



which shows that x ( w ) is a linear electric susceptibility 30|, |33| . In this way, we can easily 
define the nonlinear susceptibilities of the medium via two, three and n-point correlation 
functions between electromagnetic field and the polarization field These correlation 

functions satisfy the fluctuation— dissipation theorem 
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In the next sections we 



apply the path-integral formalism to calculate the Green's tensors and Casimir energy in the 
presence of a nonlinear dielectric medium. To simplify the problem we work in Fourier space 
and separate the Green's tensor to transverse and longitudinal parts. The longitudinal part 
of the Green's function does not lead to any force since this par t of the Green's function 
is local and the local fields do not lead to any Casimir force [29|, so we only consider the 
transverse part and for simplicity drop the superscript 



IV. NONLINEAR GREEN TENSORS 



Now let us consider the contributions induced from a nonlinear medium. To elucidate the 
method, we start with the term and generalize it to include the higher order terms like 



v 



(3) „(4) 



z/ \ etc. In analogue to the linear medium we define the nonlinear Green's function of 



the electromagnetic field up to the first order of nonlinearity as 



.(2) 



G EM {x-y) 



5 2 Z[J] 



(28) 



SJem(x)5Jem(u) j ~ ' 

The nonlinear Green's function (I28p . after some lengthy but straightforward calculations, 
can be obtained in reciprocal space as follows 



«(2) 

G EM (k,u)) 



= G EM (k, u) + G EM (k, u})A^ L (u})G EM (k, u) 

+ ^(MjASH^tk^iAj'H^fk,,) + ... 



tf- W *[e( W ) + A«( W )] 



where 



A (2) 



" LKU) ~ 2^ 



(29) 



A2 



10^ 



UJ 
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«0+' 
(30) 



is the first order correction to the Green's function of the electromagnetic field due to the 
nonlinearity of the medium. In a similar way, we find for 

«(3) 1 

Gem^i^) = — or . : .(2), : : < 3 ) , N1 ( 31 ) 

where 

POO POO POO POO POO 

^ n Jo Jo Jo Jo Jo 

1 1 1 

X u\ -(u- u' 2 f - *0+ X u 2 2 - u'i - i0+ X col - u'i - *0+ (32) 

following these procedure we find the nonlinear Green's function of the electromagnetic field 
up to the n'th order of the nonlinearity as 

1 

GEM(k,w) = S^wT4hT4(^^ (33) 

where 

^ roc poo poo poo poo 

An{{uj) = (— )" 1 / dwi / du;^ / du 2 --- du' n / dw n [z/ (n) (u;i, u; 2 , ■ ■ • ,w n )] 2 
^ Jo Jo Jo Jo Jo 

1 1 1 

; ; r= X — ^ ^ X • • ■ X 



<) 2 -*0+ w 2 -u; 2 2 -*0+ w2-w«-«0+- 

(34) 

In a similar way we can obtain correlation functions among different fields. For example, 



3-points Green's function 35] up to the first order of nonlinearity is defined by 



r {2) (i ' \- S " Z ^ i 

Csem,emA* x ' x * x '^) - 6J EM (x)5J EM ( Xl )5Mx>y j = - 

= -U)iu 2 x { > (u 1 ,u 2 )G EM (x-yi',u 1 )G EM (x 1 -x',U2) 

(35) 

where 



jo Jo 2! u'l - u\ - 10+ u'l-(4-i0+' 

(36) 

is the second order susceptibility. Also we note that in deriving Eq. (l35j) only leading terms 
have been kept and very small terms have been ignored. In order to make contact with 
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standard notation, let us recall the definition of the second order of the polarization within 
the framework of response theory 



P (2) (t) = / / dtdt' X i2) (t-t',t-t")E(t')E(t') +P% ) (t), (37) 

J — oo J — oo 

where ~P N , is the noise operator up to the second order of nonlinearity. It is easy to show 
that Eq. fl36p in time-domain can be written as 

X {2 \hM) = r d^\u Jl ^M 1 \^> {1 \^)^^ x (38) 



which is consistent with the results have been reported in 24j. Following these procedure 
we find the susceptibility up to the n'th order as follows 



OO /"OO 



Jo Jo 

1 1 

x — 7j x • ■ ■ x — k -, (39) 

which in time-domain can be expressed as 

/*oo /»oo 

Jo Jo 

sinwiti sinu n t n , . 

x x ■ • ■ x . (40) 

Now if we are given definite nth order susceptibility of the medium then we can inverse the 
relations ( HOI) and using Eq. (l26|) the corresponding coupling function v^ 1 ' can be found as 



where 



(u h u 2 , ■ ■ ■ , u n ) = w , (41) 

(2^/VMx (1 K^i)]Mx (1) M] • ■ ■ imix^M] 

poo poo 

Imx {n) {oJi , ■ ■ ■ , u n ) = / dti--- dt n x {n) (h , ■ ■ ■ , t n ) sin utfi • • • sin u n t n 

Jo Jo 



UiU 2 •••Un 



du[--- du' n x {n) (u' u --- ,u' n ) 
o Jo 



(2?r) n 

1 1 
x^ 7j -X---X— ^ -. (42) 

From Eq. (HIT) it is clear that the nth order coupling function is an odd function in frequencies. 
Therefore, the nth order coupling function will be an even(odd) function if n is even(odd). 
Using this symmetry, Eq. (1341) can be rewritten as 



A( n), x 2vrz [°° f°° (lm[ X ^(u-u 2 w n ,w 2 ,... ,uy)]) 

AWr.iu) = 7- : / duo-- - / aw, 



■ Vi '°^) n Jo Jo ^ n Im[x^{u-u 2 u n )]Im[ X ^{uJ 2 )}---Im[x^{u n )\ 

(43) 
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V. CALCULATING THE CASIMIR FORCE 



A. General formalism 

In this section we calculate the Casimir force for the simplest configuration consisting of two 
prefect conducting plates separated by a nonlinear medium of width h. For this purpose we 
consider the electromagnetic field as transverse magnetic(TM) and transverse electric(TE) 
modes. These modes are represented by scaler fields which satisfy the Dirichlet 

<p(X a ) = (44) 

or Neumann 

8 nV {X a ) = (45) 

boundary conditions where X a , (a = 1, 2) is an arbitrary point on the conducting plates. To 
obtain the partition function from the Lagrangian we apply the Wick's rotation, (t — > it) 
and change the signature of space-time from Minkowski to Euclidean. The Dirichlet or 
Neumann boundary conditions can be represented by auxiliary fields i() a (X a ) as 



%(*«)) = / V[MX a )]e^ (46) 

and 

5{d nV {X a )) = J Vty a {X a )yS d ^ x ^ x «\ (47) 
Using Eqs. (146(47]) the partition function can be written as 

Z D = Z~ l I V[<p] l[V[MX a )})e s ^ (48) 

^ 0=1 

r 2 

Z N = Z' 1 / V[<p) \{V[i, a {X a )})e s ^ (49) 

J a=l 

where Sd(ip) and Sn((p) are defined respectively by 



2 

S D [<p] = J d^x{C(ip(x)) + ip(x)J2 j d^X5(X-X a )Mx)}. 



(50) 



and 



2 

S N [tp] = I S n+ ^x{C{ip{x)) + ip{x) [ d {n) X5(X - X a )d n ^ a {x)}. (51) 
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By comparing the Eqs. (I50| |5TT) and (Q, we can rewrite Eqs. ( |50l [51]) as 

2 2 

= / f[V[i; a (x)]Z(J2 [ d»X8{X - X a )1> a (X)) (52) 

a=l a=l J 



and 

2 2 



Ziv= / n%(ar)]z£ / d n X5(X-X a )d n i(j a (X)) (53) 

a=l a=l ^ 

2 2 

where Z(£ / d n X a S(X - X a )ip a (X)) and Z(J2 J d n XJ{X - X a )d n ip a (X)) are the gener- 

a=l a=l 

ating functional of interacting fields defined in Eq. (1181) with imaginary time. From Eqs. fl52l 



531) and (TT9]) the respective partition functions can be written as 

Z D = f n^tiW]^** 1 , (54) 

J a=l 



and 



Z N = J l[V[MXa)]e lSN ^\ (55) 



a=l 

where 

iS D (ipi, tp 2 ) = iJd Q Jd, (56) 

and 

iS N {il) u i/j 2 ) = iJ N Q Jn, (57) 
and Jd(X) and Jn(X) are respectively defined by 

J D (X) = J d n X 5(X - X a ) (58) 

J N (X) = J d n X5(X - X a ) d n MX). (59) 

We will define the function Q shortly. The partition functions defined by fl54|) and fl55l) are 
calculated straightforwardly 

1 



'D 



V / d^tT T ^(x^70 
and 



^det T N (x,y, h) 
where 



(60) 



(61) 



r D (x,y,h) = I \ (62) 

G(x-y,h) Q(x-y,0) 
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and 

r ( h\ r-W^-y.O) -d 2 Mx-y,h) 

T N (x, y,h) = I (63) 

l -d 2 z g(x-y,h) -d 2 z g(x-y,0) j 

where Q is the Greens function of the fields after a Wick's rotation. Now in order to calculate 

the Casimir force we define the effective action as 

S eff = -i\nZ D [h] (64) 

from which the Casimir force can be obtained easily as 

F= dS j{ {h) . (65) 
oh 

It is easy to show that the Dirichelet and Neumann boundary condition are formally the 
same and when the plates are complete conductors lead to the same result. To save brevity 
in what follows we only calculate Casimir force for Dirichelet boundary conditions and its 
generalization to Neumann boundary conditions is straightforward. 



B. Linear and nonlinear effects of the medium (T=0) 

To calculate the Casimir energy in presence of a nonlinear medium boundary conditions are 
imposed on the electromagnetic field and the Y tensor can be written as 

„ , M r0EM,EM(x-y,O) Gem,em{x -y,h) 

LEM,EM{X,y,h) = [ \. (66) 

Gem,em(x - y,h) Gem,em(x - y,0) 

Now to obtain the Casimir force we proceed in Fourier-space since the T tensor is diagonal 
in this space. The Fourier transformation of Qem,em(x — y,h) is 



GEM,EM(p,q,h) = J dx J dye ipx+iq - y G EM ,EM(x -y,h) 

-h^[^ P0 )+^\(t Po )+---+A^l{i P0 )]pl+pl+pl 



(27rYS(p + q), 



2yJ[e(%Po) + A® ( tP0 ) + • • • + A&> (ipo)]p3 + V \ + vl 

(67) 

where p — (p , p) and p is a vector parallel to the conductors, p is the temporal component 
of p. Thus, in this case, the Casimir force is 

f - f d3p r Q{p) i r68^ 
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where 

G(p) = ^P 2 + SS^) + ^nlM + ■■■ + A^( IP0 )]. (69) 

It is easy to show that the Casimir force in the presence of a nonlinear medium ( ffTI) is similar 
to the Casimir force in the presence of a linear one, the only difference is in the definition 
of Q(p). In the absence of the nonlinearity i.e., A^ l {oj) = ■ ■ ■ = A^ L (u) = 0, the original 
Casimir force between two plates in the presence of a linear medium is recovered. 



C. Finite temperature 



Our considerations so far have been restricted to zero temperature. In fact the temperature 



A |39|-|40| where 



corrections to the Casimir force turned out to be negligible in experiments 
the measurements were performed in the separation range h < 1mm. But, at h > 1mm, as 
in Q], the temperature corrections make larger contributions to the zero-temperature force 
between perfect conductors. The generalization of the formalism to this case is straightfor- 
ward. The inclusion of temperature may be done in the usual manner 3J, 4l|]-43|]. The 
finite temperature expression, can be found by replacing the frequency integral by a sum 
over Matsubara frequencies according to the rule 



^f(iO k B Tj2Mi), £i = 27rk B Tl/h 

% 1=0 



(70) 



where T and k B are the temperature and Boltzmann constants and the prime over the 
summation means the zeroth term should be given half weight as is conventional. Therefore 
the Casimir force at finite temperature can be expressed as 



k B T 



oo ' „ 
1=0 J 



Q{p) 



2Q(p)h 



(71) 



VI. CONCLUSION 



Based on a canonical approach and using path-integral techniques, electromagnetic field in 
a nonlinear dielectric is quantized and Casimir force, in the presence of a nonlinear medium, 
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at 



finite temperature, is calculated. 
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